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Spacecraft  Potential  Calculations  — A Model 

I.  INTRODUCTION 

The  calculation  of  the  potential  on  a satellite  immersed  in  a plasma  is,  at  best, 

a difficult  problem.  In  particular,  the  accurate  prediction  of  the  potential  on  a 

spacecraft  in  the  space  environment  requires  ultimately  the  simultaneous  calculation 

of  the  paths  of  all  charged  particles  in  the  vicinity  of  the  spacecraft,  in  principle 

this  is  possible,  in  practice  it  is  not  feasible  and  a variety  of  techniques  have  been 

1 9 3 4 

developed  to  simplify  the  problem.  ’ ’ Although  capable  of  a fairly  accurate 

treatment  both  in  time  and  space,  these  models  are  limited  in  usefulness  as  they 
require  large  amounts  of  computer  time  or  do  not  include  all  of  the  various  current 
sources  necessary  to  simulate  the  charging  phenomenon.  This  report  will  present 
an  approximate  solution  to  the  problem  that  yields  spacecraft  potentials  by  making 
assumptions  that  are  equivalent  to  the  "thick  sheath"  probe  solution  for  a sphere. 
Though  similar  models  have  been  developed,  3 ’ ^ ' ' ’ 'none  have  included  the  actual 
measured  spectra  in  their  calculations.  The  model  and  FORTRAN  listings  (see 
Appendix  A and  B)  to  be  described  in  this  report  use  actual  ATS-5  and  ATS-6  data 
and  are  adapted  from  methods  originally  developed  by  Whipple  and  Del  orest. 

(See  References  1 and  9.  ) The  model,  limited  somewhat  in  its  range  of  applicability 


J 


(Received  for  publication  4 May  1978) 

(Because  of  the  large  number  of  references  cited  above,  they  will  not  be  listed  here. 
See  Reference  Page  25  for  References  1 through  9.  ) 


to  potentials  ol  -10,000  \ and  plasma  temperatures  between  50  e\  and  '10  ki-\  , 
results  in  significant  savings  in  computer  time  over  more  complicated  models. 

In  the  first  part  of  this  report  the  model  will  be  formulated.  The  individual 
current  sources  are  presented  and  approximations  as  a function  of  satellite  potential 
leveloped.  The  model  is  calibrated  using  actual  ATS -5  observations  of  the  potential. 
Two  examples  will  be  discussed  in  which  the  model  is  used  to  study  the  effects  of 
the  time- varying  geosynchronous  environment  and  of  a time-varying  photoelectron 
flux  on  spacecraft  potential.  In  Appendix  A.  we  list  the  FORTRAN  programs,  md. 
for  the  reader  new  to  the  field,  a brief  synopsis  of  she  «th  theorv  is  giv  en  in  Appen- 
dix B. 


2.  MODt.l.  KOKMl  LVTIO.N 


In  solving  the  spacecraft  charging  problem  we  are  concerned  with  finding  the 
spacecraft  potential  o such  that 


J - <J.+ J - J _ , - J . a-J.l=0 

e 1 se  si  BSe  pn 


(1> 


where 


Incident  electron  current. 

Incident  ion  current. 

Secondary  emitted  electron  current  due  to  J 


si 


= Secondary  emitted  electron  current  due  to  J 


V 


BSe 


= Back  scattered  electron  current  due  to  J , 


Photoelectron  emission  « independent  of  J or  J{(. 


Given  the  incident  ion  and  electron  particle  spectra,  the  currents  J , -Ij,  J . J ., 
and  J|^e  are  found  and  adjusted  by  varying  o on  the  spacecraft  until  Fq.  ill  holds. 
As  an  example  of  the  calculation  of  the  currents,  consider  the  following.  The 


number  flux,  XF>  , of  electrons  at  0 spacecraft  potential  is,  for  an  isotropic 


Maxwellian  particle  distribution, 

, n i2kT  1 

2 iv  = _£  ^ 

2z  1 » m ' 


\>>  = I iV*  f V 

eo  e 


1 2 


•2» 


where 


1 


.1  2 


•• • ( “:0 


in  \ 

y 

2 k T 


'•t  = mass  of  electron. 


\ = velocity, 

n(  = number  density  of  electrons, 
k 1' ^ = thermal  energy  of  electrons. 


<3> 


(See  Reference  10  :ur  letails.  > 


\l> 


eo 

qW 


This  is  equal  to  J , the  ambient  current, 
eo 


(4> 


l or  no  potential,  is  ti  . ind  Maxwellian  listribution,  Eq.  (2)  is  exact.  On 
-louever,  a probe  or  satellite  is  placed  in  the  plasma,  the  distribution  function  must 
be  altered. 

To  demonstrate  the  effects  of  ..  potential,  assume  by  I.iouville’s  theorem 
flV)=PlV’l.  (5) 


If  a is  the  charge  on  a particle 
o j 2 

! 2 m \ = 7;  m * v ' ) - qo. 

e 2 e 

r a Maxwellian  distribution, 
-ao/kT 

f'(\’)=e  ef<\> 


where 

1/2 

( j ' \ for  qo  <0  , 


O'-V  s*  for  qo  > 0 . 


For  - negative  potenti  il  <qo>0>, 

-qo/kT 

NFe(©)>=  e ' M . - (6a) 


10.  Garr--tt,  li.  I:.  • 1 :*T 7 1 V-ieling  of  the  Geos%nchronous  Orbit  I’I~rr.»  Etv.  ironm'‘nt. 
Part  i,  \i  (.!  - ; -77-0288,  ’ 


J 

j 
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For  a positive  potential  (q0<O)  the  result  is  more  complicated 


<NF  (0)  > = I IV'H'  (V”)  V'  dV1 
e 0 

00  -q0/kT 


(^) 


1/2  e 


e f(V’)  V'3  dV' 


= <NF>eo  (1+1#) 
e 


or,  converting  to  current 


je<0)  = 


-q0/kT 

Jeo  e 

Jeo(1+T3r-)  q*<0- 

e 


At  this  point  it  should  be  noted  that  Eq.  (5)  is  not  entirely  correct.  In  fact,  Eq.  (7) 

indicates  that  as  T-*0,  for  q0  <0,  J (0 ) -> <*>,  The  cause  of  this  error  is  that  in  Eq.  (5) 

the  detailed  distribution  function  F (V  , V , V ),  for  f (V),  should  have  been  used 

x y z 

(see  Appendix  B).  Along  a particle  trajectory,  Liouville's  theorem  still  holds, 

F(X,  Y,  Z,  V , V , V ) = F'  (X’,'  Y’,  Z\  V’  (0),  V’  (0),  V1  (0))  , (5’) 

x y z x y z 

where  F1  is  the  detailed  distribution  function  at  the  surface  of  the  satellite.  If  F is 
isotropic,  F’  is  not  in  general  isotropic  except  as  an  approximation.  Further,  the 
integration  of  F\  as  in  Eq.  (6),  is  extremely  difficult  for  now  the  velocities  V , Y^, 
and  are  related  to  V'x>  V'  , and  V'z  by  complex  formulas.  In  our  model,  recog- 
nizing that  Eq.  (5)  is  inherently  inaccurate,  we  will  nonetheless  assume  it  leads  to 
approximately  correct  results  (that  is,  it  yields  the  results  for  a thick  sheath). 
Eventually,  we  will  test  this  assumption  bv  comparing  predictions  of  the  model  with 
actual  data. 

Assuming  Eq.  (5)  to  hold,  this  implies  that  the  distribution  function  measured 
at  energy  E at  0 potential  is  shifted  in  energy  by  ± |qo!  depending  on  the  sign  of  the 

charge  and  the  satellite  potential.  Figure  1 illustrates  this  behavior  for  a Maxwellian 

2 

distribution  as  measured  at  the  spacecraft,.  Equating  E = 1/2  m V gives 


f(E)  = (2TTr' 


n . -E/kT 

Tm'  e 


10 
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so 


I.n  (f  (K))  - a 1 • b 

where 

a - -1/kT 

3/2 

h = Ln  I <y'7"k  ' * s'/g*  1 ’ 

hence  the  tinear  behavior  of  l.n  (ft  vs  E. 


Figure  1.  Particle  Distribution 
Function  f (assumed  to  be  Maxwellian) 
vs  Energy  E For  Various  Spacecraft 
Potentials 


For  a positive  (negative)  spacecraft  potential  and  for  ions  (electrons),  the  nar- 
ticle  distribution  function  is  shifted  to  the  left  by  energy  qo  as  shown  in  Figure  1. 
Particles  with  energv  less  than  qo  do  not  reach  the  spacecraft,  Similarly,  for  a 
negative  (positive)  spacecraft  potential  and  for  ions  (electrons),  the  particle 
listribution  function  is  shifted  to  the  right  by  energv  qo  . Since  no  particles  exist 
with  F <0  nitiallv,  the  region  0 to  qo  is  void  of  particles. 

(liven  the  nroper  integration  limits,  tin  nrobler  of  calculating  .1  , T„,  .!  . .1 

e 1 so  sl 

and  dj^e  is  reduced  to  replacing  f (I!)  bv  exp  (^^>  f (E).  In  the  following  sections 

the  integrations  leading  to  .1  , ,T,  .1  , .1  and  .1,,  will  be  performed  using  this 

e t se  SI  BSe 

replacement.  Therefore,  within  the  accuracy  issumption  >f  Eq.  (5)  (or, 

equivalently,  the  thick  sheath  approximation),  these  calculations  reflect  the  actual 
currents  to  the  spacecraft. 

Once  the  currents  to  and  from  the  spacecraft  as  a function  of  o are  known,  an 
iterative  procedure  in  which  the  spacecraft  potential  is  varied  is  employed  to  find 
0 such  that  Fq.  (1)  is  satisfied.  A FORTRAN  program  for  accomplishing  this  is 


1 1 


given  in  Appendix  i along  with  various  programs  for  calculating  the  necessary 
currents.  Input  is  required  in  the  form  of  the  differential  energy  flux  spectrum 
(obtained  for  an  electrostatic  analyzer  by  multiplying  the  count  rate  by  a constant) 
versus  energy  (62  steps  plus  one  for  background  corrections  are  employed  in  this 
version).  Details  of  the  actual  programs,  data  filtering  steps,  and  extrapolation 
are  covered  in  the  appendices. 


!.  H.M.TKON  \M>  ION  INCIDENT  C.LKRENT 


In  h'q.  (61,  the  formulas  necessary  to  calculate  the  incident  electron  and  ion 
currents  are  given  in  terms  of  the  distribution  function.  The  ATS-5  satellite,  how- 
ever, employs  electrostatic  analyzers.  Instead  of  the  distribution  function  as  a 
function  of  energy,  these  detectors  return  the  Inferential  energy  flux,  S 1 ^ * . 
Thus,  it  is  convenient  to  express  the  integrals  in  terms  or  F and  1 * ratherthan 

f (V)  and  V. 

■ r r . d (EF) . 

fne  conversion  from  f to  is  given  bv 


d (EF)  1/2  m 


Q /?  A / T-*»  y-c  t 

where  the  conversion  factors  (K  and  K.)  are,  if  f is  given  in  sec  /km,  L.  in 
,2  e 1 It 

ergs  cm  -sec-sr-e\  , and  F in  e\  : 


K - 0.  1617  for  electrons 
e 


Kj  = 5.45  X 10  for  ions 


This  gives 


<NF  > K f d <EF)  : :*L 
e e tTE  e E 


^ ^ r d (EF)  I dE 

<nfi>  = Ki  \ — 


where  the  results  are  left  in  terms  of  the  number  flux  rather  than  actual  current 
tensity  .1  (remember  only  a multiplication  bv  qs  is  required).  The  values  of 
* -jjL'are  shifted  in  notential  before  the  integration  by  utilization  of  Eqs.  (5)  and  (0). 


In  the  actual  case  an  interpolation  is  necessary  as  the  desired  value  of  d ILL'  for 
the  shifted  spectrum  usually  does  not  correspond  to  an  observed  value  of  ! ILL! 
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As  the  ATS -5  data  are  lor  liM-rete  ener  v hau  ls  for  the  r *ng»-  5t  • 51  • e\  , , 

integral  becomt  bi  mi  IE)  over  1 

these  approximations  are  a {equate  tor  a range  of  !><J  e\  to  to  ke\  in  t«  r.  per  at  u re 
and  '00  V to  -10,000  V in  potential. 


t.  Si  COM)  Alt  A I MISSION  ( l KUI  VI 


Electrons  and  ions  striking  the  satellio-  surfa  e are  either  scattered  off  the 
surface  or  cause  the  emission  of  low  end-  y,  sec  » i ,•••  . i.  itrons.  Secoi  hrv 
emission  is  usually  given  in  ten-  s >f  the  in  i iertt  differential  -urrent  dens,  v, 
dJ.  . This  current  density  is  related  to  I (Kl'l  bv 

7rr~ 


d.T. 

i 


qJ7  d (EF> 

TT  TTTT  • 


<11) 


(the  qJT  factor  will  be  dropped  in  future  liscussions  t. 

For  secondary  emission,  the  amount  of  secondary  current  et  tted  I - 
incident  flux  is  expressed  as  a ratio  5 (El, 

dJ  (E) 

6(E)  ,l.r"TP)  <>2) 

l 

where 

6 (E)  - the  yield  function, 
d.J 

- - j p-  the  differential  secondary  current  density, 

d.E 

- differential  incident  current  density  . 

The  energy  spectrum,  g(E'),  of  the  secondary  electrons  is  approximately  independent 
of  the  incident  particle  energy  and,  for  aluminur  , is  given  in  Figure  2.  * Multi- 
plying g(E’)  by  the  secondary  current  density  .1  (E)  gives  the  differential  current 
density  of  secondary  electrons  due  to  particles  of  energy  E as 


dJ  (E1,  E) 

— rrr 


g (E')  J (E)  . 


( 12) 


1.2 


SECONDARY  ELECTRON  ENERGY  (eV) 


Figure  2.  Energy  Spectrum  For  Secondary  Electrons 


From  Eqs.  (11)  and  (12) 


d J (E) 
s 


dJ.  (E) 


TTT 


= 5(E) 


6(E)  d(EF) 
E~  ~Trr_ 


This  implies  that  the  total  current  density  is  given  by. 


or  or 


'ts 


J dE'  J dE 


d J (EV,  E) 


0 0 


w 


0O  oo 


r.  d (dJ(E't  E) ) 


dE'  i.  dE  XT.  ~TT 


0 0 


OC  OO 


dE’ 


0 


r dE  d(g(E')  Js  (E)) 
0 dE 


00  «>  dJ  (E)  dE 

J g (E')  dE’  / 

‘0  0 


7TFT 


00  00  d.l.(E)  dE 

Jo  ^E')dE'in  6(E)-W- 


I g(E')  dE'  J 6(E)  [^£^1  JJf-  . 

‘0  n L nr.  j t 


0 


(14) 


(15) 


The  energy  spectrum  g(F')  peaks  sharply  at  ~2  eV  and  thus  for  ions  and 
electrons  is  ~ 0 if  0 is  greater  than  • 2 \'  (all  the  secondary  electrons  are  attracted' 
and  1 if  0 is  negative.  To  the  order  of  accuracy  of  these  calculations 


j (I-  ' 11  ' will  be  assumed  - 0 for  o * 0 and  1 for  0^0.  ! rort  data  presented 

* * 1 

by  Whipple  0 ) is  approximated  for  electrons  :.»•  mu-tin:  <>n  aluminun  bv 
6 (F)  ; 0.  018  ' fl - A 1 


0.786  tog]0  }•:  -o.  a'is 


- 0.  45  (log  E)  • 2.  1 
0-34e-K/28800 

The  aetual  curve  is  shown  in  Figure  8. 


1 •'  20  e\ 


20  e\  < F < 100  c\ 


(16) 


0.  18  sin  I (log,-  E - 2.)]  ■ u.75  100  e\  E * 10  e\ 


10,?  eV  < F ^ 2X  104  e\ 


2X  10  eV  < 1.  . 


Figure  2.  Secondary  Electron  Yield  for  Electron  Impact 
>n  Aluminum 


I he  vield  function  b (K)  is  not  well  known  for  ions  impacting  on  aluminum  but 
,i  it  to  the  lata  presented  bv  Whipple  (see  Figure  4)  indicates  that  (or  H and 
for  the  energy  range  of  interest 

<5  (F)  * 0.  086  E < 700  eY 

’ ‘ "*P  (ttStot)  700  PV  ? 5 5°5  “V  • (V 

These  values  should  be  used  with  caution. 

Note:  We  are  onlv  interested  in  the  energy  range  0 to  85*  keY  in  this  study.  The 
actual  rnrve  turns  over  at  100  ke\  . 
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Figure  4.  Secondary  Flectron  Yield  for  Ion  (II  ) 
Impact  on  Aluminum 


The  yield  function  6(E)  for  the  incident  ions  an'1  electrons  is  substituted  into 
Kq.  (15)  and  integrated  using  the  appropriate  shifted  values  for  1 1 to  obtain  the 

secondary  emission  currents  due  to  the  electrons  and  ions.  Typically,  the  second- 
ary emission  due  to  electrons  results  in  a current  of  approximately  25  to  50  percent 
that  of  the  incident  electron  current  while  the  secondary  emission  due  to  ions  re- 
sults in  an  electron  current  approximately  2 to  3 times  that  of  the  incident  ions. 

As  will  be  discussed  shortly,  it  was  necessary  to  adjust  these  values  to  obtain 
accurate  estimates  of  the  potential. 


5.  HACkSC  VITKREI)  KI.KC.TRON  CURRENT 


Some  incident  electrons  are  reflected  giving  rise  to  the  backscattered  electron 
current  (the  backscattered  current  due  to  inns  is  very  small  and  ignored).  Although 
-iri  pie  theories  ol  collisional  scattering  such  as  those  of  Everhart  are  useful  in 
predicting  the  net  current,  experin  ental  curves  for  backscattered  emission  from 

12  <i 

the  intlt;  are  also  available  * St  e rnglass ),  “ DcEorest  lias  adapted  these  results 
to  the  ATS-5  lata.  Ilis  levelopn  ent,  which  v.  ill  be  uresente  I here,  is  similar  to 
that  given  'or  secondary  electrons. 


11.  1 r.  K«  I960)  f pie  theor’  t < t i t tiot 

solids,  rna  f Vppl.  PI  v ■.  31  . >.  8 1;  I 183. 

1 2 . SI  ...  1 1<  : t 

Phvs.  Rev.  '‘5:345, 

- — ■ ■ — 
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The  basic  equation  for  baekseattering  is 


« <1 .1.  (i:> 

I B(E\  K)  -jJ-  dj; 
I • ltTr 


F'  - energy  of  backscattered  particles  ( 1 ; ’ ? 1), 

I ■ energy  of  incident  particles, 

B(F',  F)  = percentage  of  electrons  scattered  at  a given 
energy  K'  as  a result  of  an  incident  electron 
at  energy  1 , 

• I.(l'i)  = incident  current  (electrons  in  this  case). 


From  Sternglass 


B(K',  1 :> 


fid.’  / i > 


O is  given  as  a function  of  K I'M  (Deforest,  private  communications)  in 
Figure  5.  Numerically,  for  aluminun 


c;  (K ) = -0.  025  Is  • 0.  925 


K " 0.  8 . 


-0.725  K + 0.7  65  0,6  fK50J 


■ 0.  05  K t 0,30  O.isKf  0,6, 


0.  5 K +■  0.  12 


K ‘ 0.4 


(Note:  G (K)  ~ 2y  (1-K)). 


For  aluminum,  when  B (E1,  E)  is  substituted  in  Fq.  (181,  a ratio  of  '25  percent  for 
the  backscattered  current  to  incident  current  is  obtained  in  agreement  with  other 
estimates. 


13,  Grard,  R.J.L.,  Knott,  K. , and  Pedersen,  \,  (1973)  The  influe nci  photo- 

electron  and  secondary  electron  emission  on  electric  field  measurements 
m the  magnetosphere  and  solar  wind,  Photon  and  Particle  Interactions  with 
Surfaces  in  Space,  R. Grard  (FI.)  I).  Reidel  Publishing  Co.  , 

I lordrecht,  Holland,  pp  163-100. 
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Figure  5.  Graph  of  11  (E',  F)‘  E As 
a Function  of  (E'/E)  Where 
B(E',  E)  Is  Approximately  the 
Percentage  of  Electrons  Scattered 
at  a Given  Energy  E'  As  a Result 
of  an  Incident  Energy  E 


Continuing,  Eq.  (18)  becomes 


• c m /*[f  o(f)]qpi.N4- 


(Note:  — p—  md  — = 0.  12  for  ATS -5) 

Substituting  the  proper  values  of  (^jji  * p in  I q,  (21)  and  performing  the 
integration  gives  the  total  current  due  to  backscatterod  electrons.  Enlike  the 
secondary  electron  current,  though,  the  backseattered  current  is  a gradual  function 
of  spacecraft  potential.  For  positive  potentials,  the  0 integral  limit  in  Eq.  (21) 
should  be  replaced  by  qt>  , the  energy  shift  lue  to  the  spacecraft  potential. 


().  PIIOTOKU  a ilON  EMISSION 


Light,  particularly  in  the  ultraviolet,  pacecrafl  ts<  • ••  is  - 

sion  of  photoelectrons.  Although  the  haractt  risti  tti  rticl<  nd  t • 

processes  involved  are  well  known,  the  actual  photoelectron  ett  ission  fror  a spaoe- 
craft  i poorly  known,  rh<  i et  1 th<  tee- 

craft  surface  arid  the  lack  of  precise  knowledge  of  ’tie  solar  -oe  trun  and  its  inter- 


iction  with  various  materials.  Gr art  . t ,1 


ombined  the  sol  r 
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SIX'1  trill!  with  the  emission  characteristic-  ,.r  various  substances  to  jive  the  photo 
1 ' 1 " irrent  a tioi  enei  . tel  est  nicatioi 

algnrit  >xi  t theirri  Its.  It  ves  the  photoelectrot 

current  is  i function  of  positiv,  spa, •♦•••raft  potential  (for  negative  potential  all  of 
the  photoelectron  ‘ u t •■< ■ t it  leave-  the  -pare  raft).  It 


i~n.  ; 


(.1  > 0 . 


(22  I 


< 1 n 1 ible  !.  i rot  table  it  ear  that  est  tes  o 

L0  range  ov,  ...  . nitude.  Fortunately.  ■ ■ 

that  we  will  be  concerned  wit  nvo  • hadowed  irfac es  ■ • , j r i 

no 

I he  charging  model  will,  however,  in  onjunctinn  with  actual  : ita,  he  ,-c  i jn  i 
later  -ection  to  cstic  ate  tin  \alne  r appropriate  to  \ is-i. 


Table  1.  Values  for  .1 

po 


Author 

2 

saturation  Current  A cm“) 

Grard  et  al 1 

4.  20 

(Alun’inun  oxide) 

00 

(Indium  oxide) 

n 

Me  Forest 

0.  40 

(Graphite) 

0.  H2 

(ATS -5) 

Whipple* 

!.  00 

(Hocket  -aluminum ) 

T.  COMPARISON  \\  Mil  D ATA 

The  ultimate  check  of  any  model  is  how  well  the  predicted  results  compare  with 
easuret  ents.  The  basic  set  of  lata  will  he  spe  ttra  fr  > \TS-5  ot  erio 

immediately  before  and  after  entry  into  the  earth’s  shadow  and  for  periods  r e li- 
ately  before  and  after  exit  Iron,  the  earth's  shadow.  These  periods  were  seleete  i 
as  they  -an  bo  used  not  only  to  study  large  potential  variations  (on  the  order  of 

tlibi  ■ ■ ie  p »1  telectron  flux.  Table  2 lists  the  eclipses  studied 
and  the  observed  potentials  during  eclipse. 
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Table  2.  Eclipses  Studied 


Date 

UT 

Potential 

1969  22  Sep 

0629 

-3400 

07  31 

-3810 

16  Oct 

0627 

-5360 

07  11 

-3810 

1970  12  Sep 

0631 

-2420 

07  18 

-1730 

15  Sep 

0626 

-877 

0721 

-1540 

17  Sep 

0623 

-5380 

0723 

-3040 

19  Sep 

0620 

-2720 

07  24 

-1940 

26  Sep 

0724 

-2170 

17  Oct 

0630 

-1230 

0659 

-2170 

18  Oct 

0633 

-397 

0650 

-316 

19  Oct 

0640 

-396 

0648 

-558 

ATS-5  does  not  consist  of  a single  material  nor  can  it  be  said  to  be  spherical 
in  shape  (a  tacit  assumption  in  the  preceding  analysis).  However,  keeping  to  the 
spirit  of  a "simple''  charging  model,  the  satellite  was  approximated  as  an  aluminum 
sphere.  Figure  6 shows  the  results  of  these  calculations.  The  liscrepancies  be- 
tween observed  and  predicted  potentials  have  been  corrected  bv  adjusting  the  coeffi- 
cients of  the  secondary  emission  terms  to  give  a best  fit  in  a least  squares  sense 
(the  backscattered  flux  correlates  so  well  with  the  incident  electron  flux  that  deter- 
mining its  coefficient  is  somewhat  difficult  as  it  may  reflect  slight  errors  in  the 

actual  measurement  of  J ).  For  correction  factors  of  1.3  (.1  ),  0.55  (.1  ,),  and 

e se  si 

0.  4 (J  g),  the  results  have  a standard  deviation  of  * 800  volts.  It  should  be  kept  in 
mind,  though,  that  the  model  based  on  several  assumptions  and,  considering 
these,  this  agreement  is  quite  good. 
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ATS -5 


Figure  6.  Predicted  vs  Observed 
Potentials.  The  results  for 
spectra  when  the  satellite  was 
illuminated  and  when  it  was 
shadowed  are  given 


8.  MOM I.  APPLICATION 


The  model  has  been  developed  for  two  purposes.  First,  combined  with  the 
geosynchronous  plasma  model  developed  by  Garrett,  it  can  be  used  to  predict 


potentials  on  spacecraft  as  a function  of  the  geomagnetic  index  A and  LT,  the 
satellite  local  time.  It  also  is  employed  to  calculate  the  potential  on  a spacecraft 


as  it  passes  into  and  out  of  the  earth's  shadow. 


Figure  7 plots  the  potential  variations  predicted  bv  the  model  for  the  geosyn- 

chrous  simulation  model  of  Garrett  ^ for  moderate  geomagnetic  activity  (A  - 120) 

P 

and  high  geomagnetic  activity  (A  1656).  That  simulation  gives  a ' 2 Maxwellian 
fit"  to  the  plasma  such  that  the  distribution  functions  are  for  electrons 
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figure  7.  Predicted  Potential  of  a Shadowed,  fleet rically  P dated 
Surface  as  a function  of  focal  Time 


The  distribution  functions  are  converted  to  differential  energy  spectra  using  f'q.  I'd. 

••  3 pi  • . in  nserted  into  the  progran  u 1,  issui  ing  no  photoelectron  current, 
the  potential  calculated.  These  potentials  are  ihe  : .aximum  that  would  be  expected 
for  a shadowed,  electrically  isolated  surface  as  a result  of  the  ambient  environment. 

Another  use  of  the  model  is  in  the  leterr  ination  of  the  potential  as  a spacecraft 
thi  earth's  sha  low.  The  eclij  i i ibed  earliei  n isi 

r in e the  photoelectron  flux  necessary  to  give  the  observed  potential  variations  as 
NTS -5  passed  into  an  , out  of  the  earth's  shadow.  By  varying  the  at'  plitu  ie  >•’  an 


appropriate  model  of  atmospheric  attenuation  to  fit  the  observed  photoelectron  <ur- 

2 

rent,  .1  was  found.  The  resulting  value  was  determined  to  be  0.  -4  n A m“  w hich 

po  14 

is  in  agreement  with  the  values  in  Table  1.  The  reader  is  referred  to  Garrett 
for  details  of  this  procedure. 

Figure  8 shows  the  observed  and  predicted  potentials  as  ATS -5  passes  into  and 
out  of  eclipse  on  17  September  1 f *7 0 . The  results  of  the  prediction  are  adequate  and 
within  the  j 800  V error  but  their  deviations  from  the  observed  values  m as  indicate 
either  a need  to  include  sheath  effects  or  that  the  thick  sheath,  spherical  probe  ap- 
proximation is  inaccurate.  In  any  event,  this  method  of  testing,  by  comparing  the 
observed  and  predicted  potentials  as  the  photoelectron  flux  is  varied,  should  prove  tobe 
a powerful  tool  for  comparison  with  other  spacecraft  charging  models  in  the  future. 


UT(mln) 


figure  R.  Cop  narison  of  Observe  I mi!'"  1 *o 

the  \TS-5  Satellite  as  a Function  of  '!  < . I ope 

potentials  predicted  by  eclipse  spe< : rn,  lower  i 
'ials  predicted  by  illuminated  -pc  r , , m i . -tr 
average  of  the  two 

: rr<  tt,  H.  B. 

craft  Potential,  A I ... 
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A simple  model  based  on  the  work  of  Whipple  and  Def  orest  has  been  develope  1. 
The  model  was  calibrated  with  ATS -5  plasma  data.  The  model  predictions  for  the 
potential  on  a satellite  as  it  passed  into  eclipse  were  compared  with  actual  observa- 
tions. l'he  results  indicate  agreement  between  the  predicted  and  observed  values. 

The  model  is  used  in  conjunction  with  a model  simulation  of  the  geosvnchrous  en- 
vironment to  predict  spacecraft  potentials  tinder  different  geomagnetic  and  local 
time  conditions.  The  model,  after  being  calibrated,  successfully  predicts  ootentials 
with  t H00  V accuracy  over  a range  of  10,  000  V.  It  includes  relevant  current  terms 
and  is  efficient  in  comparison  wdth  other  more  complex  models,  taking  '2  sec  per 
potential  calculation.  A FORTRAN  listing  of  the  program  is  provided  in  Appendices 
A and  11.  For  reference,  a simple  theory  of  the  plasma  sheath  is  also  presented. 
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Appendix  A 

FORTRAN  Listings 

The  solution  of  Eq.  (1),  the  current  balance  equation,  involves  an  iterative 
procedure  with  three  steps.  Once  the  input  data  are  properly  filtered,  the  interation 
commences  with  a computation  of  the  currents  assuming  the  initial  potential  to  be  0. 
The  currents  are  combined  as  in  Eq.  (1)  and  the  deviation  from  0 current  calculated. 
Depending  on  whether  the  deviation  is  close  to  0 or  not,  the  process  is  either  ter- 
minated or  a new  potential  estimated.  The  spectra  are  shifted  by  the  new  potential 
and  the  process  repeats.  (See  flow  chart  at  end  of  this  Appendix.  ) 

As  the  determination  of  the  energy,  E,  and  of  the  corresponding  differential 
energy  flux,  1 j - , for  the  electrons  and  ions  is  satellite  dependent,  we  will  not 
discuss  this  aspect.  The  program  as  currently  arranged  accepts  an  energy  spec- 
trum having  N steps.  The  first  step  is  assumed  to  be  0 energy  and  contain  the 
background  counting  rates.  The  subroutine  listings  using  this  as  input  follow 
Appendix  B. 

The  'lata  from  a real  satellite  are  usually  not  continuous.  Subroutine  1TL  is 
designed  to  take  the  dat  ■,  interpol  te  missing  points,  drop  values  below  background 
or  below  the  satellite  potenti  >1  < p r t i < les  h iving  energy  less  than  the  spacecraft 
potential  are  trapped  and  do  not  contribute  to  the  satellite  current)  and,  if  desired, 
estimate  values  outsi  !e  the  energv  r mge  ■ >f  the  instrument.  In  the  current  form  the 
program  ignores  .11  channe  ls  . inr  0 < . r neg  itive  energy  flux  until  the  first  non- 
zero channel  i-  re  i e 1 , , Vlt.  r this  channel,  ill  mis-ing  values  are  interpolated 
using  a 1 ■ . the  last  ro,  it  will  be  left  zero.  The 


I 


program  exports  \ values  of  tin  ■ lilTcrent  i.d  energy  flux  in  on  » i 1 •••loot  runs  I 
of  e\  1 < ■ • ■ ■ ' ■ ■■  ' ■ : ! ’ 

tin-  measurements  wore  made  (held  in  iwnnwtil. 

tine  -1111  I take  s is  input  Cl,  C2.  and  E.  It  also  nei  the  i ’ < ■ 

• , . , ind  the  satellite  number  (i  iportant  only  for  VI  >).  I 1 ' ' ' ■ 1 

which  tile  spectra  wore  determined  (tint  is,  spacei  r ill  potential!  is  \ S.  1 he 
potential  to  which  the  spectra  ire  to  be  shifte  i is  V 0.  rhe  new  value  of  Cl  nd  C2 
are  found  by  first  determining  the  energy  in  the  unshifted  spectrum  corresponding 
to  the  new,  desired  (shifted)  energy.  SI  (or  S2),  the  shifted  value,  is  etei  ined 
by  interpolating  between  tile  \ dues  of  Cl  (or  C2)  eorr.  spoil  ding  to  tin  - in 
bracketing  the  unshifted  energy.  Using  I-iouville1  s theorem,  this  Interpol  te 
is  shifted  to  the  new  energy.  E , by  dividing  it  by  the  square  of  th<  rati  f thi 
unshifted  energy,  K , to  tile  new  energy 

si  (i;  ) - cl  (i:  i (i : / 1 : r . <Al> 

s o S () 

The  results  of  Sllll  1 are  input  to  subroutine  1)11  1 where  the  in  lr.  i hi  d i urrents 
ire  i il<  ul  ted.  The  incident  currents  due  to  the  electrons  ' I id  in  i ions  1 1 I • ire 
calculated  in  TOTAL  according  to  Eq.  (10)  tnd  return!  ■ 1)  • x(2)  (all  cur- 

rents  are  in  n/(cm"-sec-sr)  rather  ttian  n Vein-).  Till  \I.  require.-  funetion 
D(E,  I),  w here  E is  the  energy  irraj  in  1 I is  t hi  encr(  stej  i • ' ! ■ 

ound.  D is  the  r itio— . A E(l)  is  the  enerj  cl  nnel  width  for  letector. 
Mthough  function  D is  satelliti  lependenl  (as  n exai  pl<  , on  VI  S-5  A E E = 0.  12), 
it  can  be  generalized  if  the  number  (id  in  e .ri  Iti  is  changed. 

] e backscatteri  lux  i lculated  BACK,  rhe  total  current,  issuming  0 
or  neg  . t i v e potenti-d,  i.-  returned  dong  witu  the  b.icksi  itti-red  current  »s  a function 
of  positivi  potenti  1 (DBSEI.  Similarly,  thi  second  ry  currents  due  to  the  electrons 
(SEi  is  total)  ind  ions  (Sl'l  is  total)  are  returned  in  Ml)  and  M2)  by  subroutine 
-i.i  , The  phol  1 1 d al  )i  1*  1 1 Eq.  (22).  Thel 

scattered  current  is  ri  turned  in  Mat,  the  secondary  current  due  to  electrons  in 
Xld).  the  secondary  current  due  to  ions  in  \(4),  and  the  photoelectron  current 
in  M6  >. 


1 iu-  secondary  flux  due  to  ions  is  not  icrurate  dmvc  100  ke\. 


Appendix  B 

Sheath  Theory 

The  measurement  of  the  characteristics  of  .i  contained  plasma  in  the  laboratory 
is,  in  conception,  quite  simple.  Basically,  a probe  (a  wire  for  example)  is  placed 
in  the  plasma.  The  potential  on  the  probe  is  varied  and  tile  resulting  current 
measured.  A typical  current  vs  voltage  curve  is  shown  in  f igure  Bl.  As  in  most 
respects,  the  probe  problem  is  identical  with  that  of  a satellite  immersed  in  a 
plasma,  our  discussion  will  begin  with  a description  of  probe  theory. 

Qualitatively,  if  the  ion  and  electron  plasma  populations  can  be  Inscribe  I b\  < 
Maxwellian  particle  distribution,  f igure  Bl  can  be  readily  interpreted  as  follows. 
Referring  to  the  figure  in  region  A,  the  ion  current  is  collected  and  the  electron 
current  repelle  1 it  the  probe  potential  is  negative.  In  region  B tin-  electron  current 
is  only  weakly  repelled  by  the  probe  potential.  At  the  border  between  region  li  and 
(',  the  ion  and  electron  currents  are,  in  fact,  equal  (this  corresponds  to  the  space- 
craft charging  problem),  in  region  C the  electron  current  dominates  the  ion  cur- 
rent (Note:  as  the  electrons  are  ItiOO  times  less  massive,  for  the  same  energy  they 
move  40  times  faster  than  the  ions  so  the  electron  current  is  ~40  times  larger). 

In  region  If  almost  all  the  ion  current  is  repelled  and  the  electron  current  saturates. 
The  exact  curves  for  all  these  regions  are  dependent  on  the  temperatures  and  den- 
sities of  the  ions  and  electrons. 


' 


I igure  HI.  Probe  < urrent  (Je)  • - 
Potenti  al  0 1 or  . Spheric.il  Prcihe 
Immerse'l  in  a Neutral,  \l  .xwellian 
Plasma 


I'lus  simple  picture  neglects  one  important  problem  — the  sheath.  As  either  tin 
. >1  electr  mt  ■ in  ite  the  current,  the>  foi  n < lo  in  the  ii  inits  of  th< 

iri'l  i n drastically  liter  the  calculation  of  the  currents  shown  in  Figure  HI.  To 
>bt  un  better  description  of  the  sheath  region,  consider  the  following.  Assume 
th  t the  electrons  in  thermal  equilibrium  ibove  eh., rge  1 pi  te  (ti.  t is,  .issue  e 
the.  form  m itmosphcre"  with  gravitv  replaced  by  the  electric  field): 

I q i 0 ( X ) 

kT 

n (X)  = n e c ( HI ) 

e o 

where 

n number  dens  it  v of  electrons  at  surface  of  plate, 

o 

n = number  lens  it  v of  electrons  is  a function  of  \, 
e 

X = height  above  charged  plate  (\  0 at  plati  surf  cel, 

O potential  (assumed  positive), 

e]  = electron  charge 

kT  = thermal  energy  of  electrons. 

t irther,  issume  that  the  ions  ire  motionless  nd  nifori  1 . 1 . i . t . B\ 
Poisson's  equation 
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( B21 


-aj  = - 4 71 1 q | (n.  - n ) , 

dX  1 e 


where 


ii-  number  density  of  ions. 

1 

.Substituting  in  for  n.(n.lx)  ~ n ) and  n , for  small  [ 'j.J  ^ 
h ii  o e k 1 t 

I q 1 0 

2 kT 
— f-  = -4ff  | q | n (1  - e e ) 

<1 X ° 

-4ff|q|no  (|q|0/kT  > . (B3) 

Thus 

-(x/X  1 

0 = 0 e ° < B4  > 

o 

w here 

/ kT  . V " 

4 = I ^ — j = Debve  length, 

V"q  n<7 

0(j  = potential  at  plate  surface  . 

This  distance  (the  Debve  length)  is  the  distance  from  the  probe  over  which  the  probe 
potential  is  attenuated  by  l/e  and  is  the  distance  over  which  the  probe  potential 
disturbs  the  plasma.  This  is  usually  called  the  sheath  region. 

The  she  th  region  around  the  probe,  particularly  if  it  is  large,  c .n  significantly 
alter  the  shape  of  the  curve  shown  in  f igure  HI.  Within  it,  particles  are  influenced 
bv  the  changing  probe  electric  field  and  detailed  single  particle  orbit  calculations 
.re  often  necessary.  In  space  the  number  densities  are  quite  low  (<  1 00/cm'  ) and 
the  temperatures  often  on  the  order  of  20  ke\  making  the  Debye  length  (or,  approxi- 
mately, the  she  'th  region)  on  the  order  of  tens  of  meters.  The  plasma  is,  however, 
collisionless  — somewhat  simply ing  the  computations. 

The  effects  of  particle  orbits  within  the  sheath  can  be  seen  in  the  following. 

A ssuming  a thick  sheath  (that  is,  that  A • • r , w here  r is  the  r idius  of  tile  s itel- 

o s i s 2 

lite),  consider  a particle  having  initial  kinetic  energy  m V()  and  approaching  the 
probe  within  a distance  r.  Then,  by  conservation  of  energy  and  ingular  moment un 
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( 1 55 ) 


1 •>  1 2 

— in  \ ” -T  Hi  \ + qo  (f), 

2 o - r 


m K V m r \ 


where  the  qq  is  the  potential  energy  (0  at  infinity)  and  m is  the  mass  of  the  particle. 
Solving  for  the  impact  parameter  K 


1{“  r“  (1  - qO  (r)  1 

T 

m V 


(HU) 


(■ 


-) 


Therefore,  the  total  current  density  at  the  surface  of  a spherical  probe  of  radius  r_, 
is.  tor  a monoenergetie  particle  population,  approximately 


J =- 


1 


4n  r 


T = Jon  -J2. 
; m V 


) 


(157) 


where 


J..  = 


q n \ 


4 77  K 

= the  current  per  unit  area, 
= ambient  density  it  K. 


This  follows  as  any  particle  which  would  normally  pass  through  a volume  of 
radius  R would,  in  the  presence  of  ( potential  o,  pass  through  the  volume  of  r 
..  , ssumed  to  be  the  spacecraft  surface).  This  would  give  a current  1 = 4r  R JQ 

at\. 

, 15 

l or  i Maxwellian  distribution,  Kq.  (Hi)  becomes 


.1  = T = .1  ( 1 


Iqoi 


kT 


where  kT  has  replaced  i m V and 


1/2 


i . ,n  t v 

-*Un  ( — : ) 

o 2 it  m 


See  Kq.  ( 1 ) in  Reference  10. 


! ( hen.  I .1  . (1965)  Klectric  Probes,  Plasma  Diagnostic  Techniques.  K.  II. 
1 1 uddestone  and  S.  I..  Leonard  (Kds.  ) Academic  Press,  Inc.  , 

New  York,  New  York,  pp  113-200. 
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This  result  is  only  applicable  for  A >>  r the  thieu 

>4  o Is’  the  thick  sheath  approximation,  and  for 

a sphere  (see  Chen  for  details  of  this  development). 

i'hese  results  for  the  spherical  thick  sheath  approximation  are  equivalent  to  the 

model  presented  in  this  report.  As  described,  explicit  assumptions  about  the  nature 

“r  PlaSma’SPaCeCraft  intera‘tion  made  that  result  in  an  expression 
Itq.  (7)]  identical  with  that  for  a spherical  thick  sheath. 

Finally,  as  discussed  in  detail  in  Reference  10  the  plasma  distribution  function. 
(X.  Z,  Vx.  Vy.  Vz)  (where  X.  V.  Z are  spatial  coordinates  and  V V . V 
are  velocity  components)  is  the  conventional  basis  for  describing  a space  plasma*. 

This  function  was  simplified  by  assuming  an  isotropic,  collisionless  plasma  which 
gives: 

FIX,  V.  Z.  Vx,  Vy,  Vz)  dxdy  dz  dV'x  dVy  dV?  = f (X,  V,  Z.  V)  dx  dy  dz  (4  * V2dV) 
where 

2 2 n 1 / 2 

v=(V  + Vy2+  vz2>  . 

As  discussed  in  the  text,  this  is  an  approximation  and  must  be  used  with  care. 
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